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Abstract 



We propose an amplitude-phase representation of the dual-tree com- 
plex wavelet transform (DT-CWT) which provides an intuitive interpre- 
tation of the associated complex wavelet coefficients. The representation, 
i— i in particular, is based on the shifting action of the group of fractional 

f~H Hilbert transforms (fHT) which allow us to extend the notion of arbi- 

trary phase-shifts beyond pure sinusoids. We explicitly characterize this 
<Z3 shifting action for a particular family of Gabor-like wavelets which, in 

i ^ i effect, links the corresponding dual-tree transform with the framework of 

windowed-Fourier analysis. 

y—i We then extend these ideas to the bivariate DT-CWT based on cer- 

tain directional extensions of the fHT. In particular, we derive a signal 
representation involving the superposition of direction-selective wavelets 
affected with appropriate phase-shifts. 
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1.1 The dual-tree transform 

T T te begin by briefly reviewing the fundamentals of the dual-tree transform. 
• »-j V V The transform involves a pair of wavelet bases with a one-to-one 'quadra- 

ture' correspondence between the basis elements [1, 2]. Specifically, one con- 
siders a primary wavelets basis {tpi,k}(i,k)ez 2 of L 2 (R) generated through the 
dilation-translations of a single prototype ip(x); that is, ipi,k{x) — ^i,kip(%) 
where Ei^f(x) = 2 l l 2 f(2 l x — k) denotes the (normalized) dilation-translation 
operator corresponding to integers i and k. The highlight of the transform is 
then the construction of a secondary wavelet basis {ip^ k }(i,k)ez 2 having the cor- 
respondence k (x) = J^fipi^k(x), where denotes the Hilbert transform (HT) 
operator: 

J?f(x) -j signH/M. (1) 

The HT acts as a quadrature transform that takes cos(wo:e) into sin(woa:), and 
as an orthogonal transform on L 2 (R) in the sense that (f,J4?f) = for all 
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f(x) in this space. Though this is not at all obvious a priori, it turns out (as 
suggested by the notation) that the secondary wavelet basis can also be realized 
through the dilations-transations of the HT counterpart ip'(x) = Jt°ip(x). This, 
in fact, is possible thanks to certain fundamental invariances enjoyed by the HT 
operator. In particular, following definition (1), one can readily verify that the 
HT commutes with translations and dilations; in particular, 

J4? Ei t k — S^ fe Jf ; (2) 

and that it is unitary: 

IWIIl2 = H/IIl2 (/ € L 2 (R)). (3) 

It is then easily deduced that the functions ip^ k (x) — Ei^ip' (x) indeed constitute 
a wavelet basis of L 2 (R), and that the correspondence ip[ k (x) — Jfipi v k(x) holds 
for every integer i and k [3]. 

The application of the transform involves the simultaneous analysis of a 
signal f(x) in L 2 (R) in terms of the quadrature wavelet bases {f/'i.fe} and {tp^ k }. 
In particular, one considers the wavelet expansions 

where the expansion coefficients in (4) are specified by the dual wavelet bases 
{V'i.fc} an d {Wi k} through the projections 

ai[k} = if,i>i,k), and bi[k] = (/, r/S- )fc ). (5) 

In effect, this allows one to identify the complex wavelet coefficients Cj[fc] = 
(ai[k]+jbi[k])/2, and the associated amplitude-phase factors \ci[k)\ and arg(cj[fc]) 
(the use of the factor 1/2 will be justified shortly). As a consequence of (2) 
and (3), the dual wavelet bases can also be generated through the dilations- 
translations of two dual wavelets, ip(x) and ip (x), that form a HT pair as well. 



1.2 Multiresolution Gabor-like transforms 

A framework for constructing HT-pairs of wavelets (within Mallat's multiresolu- 
tion framework) was recently proposed based on a spectral factorization result 
for scaling functions [3]. In particular, it was shown that a multiresolution 
form of Gabor-like analysis could be achieved within the framework of the dual- 
tree transform. This was founded on the two vital observations. The first one 
was that the extended (a, r) family of B-splinc wavelets tp{x;a) (indexed by 
the approximation order a + 1) is closed with respect to the action of the HT 
operator (the corresponding discrete wavelet transform has an efficient FFT- 
based implementation) . Secondly, it was shown that the complex spline wavelet 
\$>(x) = ip(x; a) + jJf?ip(x; a) asymptotically converges to a Gabor function: 

*(a;; a) ~ <p(x) exp(juj x + £ ) (a — > +oo), 
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where <p(x) is a Gaussian window, and u>o,£,o are appropriate modulation pa- 
rameters. As a consequence, a Gabor-like transform, involving the computation 
of the sequence of projections 

f(x)»±(f(x),Ei tk *(x;a)) (i.fceZ) (6) 

with the dilates-translates of the Gabor-like wavelet ^f(x;a), could be realized 
using the dual-tree transform corresponding to the spline wavelets tp(x; a) and 
J4?ip{x; a) (sufficiently large a). These ideas were also extended for the realiza- 
tion of a direction-selective Gabor-like transform. 

1.3 Present Contribution 

In this paper, we provide a characterization of the dual-tree transform, and 
the Gabor-like transforms in particular, from the perspective of multircsolution 
windowed- Fourier analysis. In particular, we link the multircsolution wavelet- 
framework of the former with the intuitive amplitude-phase representation as- 
sociated with the latter. Complex wavelets, derived from the combination of 
non-redundant wavelet bases, provide an attractive means of encoding the rel- 
ative signal "displacements" using the phase relation between the components. 
The DT-CWT is a particular instance where the components are related through 
the HT. 

In §2, we derive a representation of the dual-tree transform using the group 
of fractional Hilbcrt transform (fHT) operators: 

ttr = cos(ttt) J - sin(Trr) JT (r e R) (7) 

(J^ is the identity operator). In particular, we are able to interpret the phase 
factors associated with the dual-tree transform in terms of the action of this 
group. However, it is the fundamental invariances (2) and (3) inherited by 
this extended family of operators that play a decisive role in establishing the 
windowed- Fourier-like representation (cf. (11)). 

The proposed windowed Fourier-like representation admits a straightfor- 
ward extension to the bivariate setting by introducing an appropriate multi- 
dimensional extensions of the HT. In particular, we arrive at a representation 
(cf. (15)) involving the superposition of the direction-selective synthesis wavelets 
affected with appropriate phase-shifts. This provides an explicit understanding 
of the phase-shift action of the fdHT operators for a the particular family of 
2D Gabor-like wavelets derived through the tensor products of ID Gabor-like 
wavelets. 
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2 DUAL-TREE GABOR WAVELET ANALY- 
SIS 



2.1 Signal representation: interpretation of the amplitude- 
phase factors 

Our objective is to derive a representation of f(x) in terms of the amplitude- 
phase factors Ci[k] — \ci[k]\e'^ i ^ . Clearly, the transformation 

f(x) ~ {ci[k]} {i k)ez2 

constitutes a overcomplete representation of f(x). In particular, given the coef- 
ficients Ci[k], there exists non-unique ways of reconstructing the input f{x). We 
consider the simplest inversion procedure involving the inversion of both the 
forward transforms, as in (4), followed by the averaging of the reconstructed 
signals. In particular, by combining the expansions in (4) and by invoking the 
dilation-translation invariance of the fHTs, we arrive at the following represen- 
tation: 

fW=l E (adk^A^ + bdkW^x)) 
(ii)ez 2 

= E l c *[ fc ]l%[fc]/7r{^fc(a;)} 

(i,k)eZ 2 

= ]T \ Ci [k\\^i,k{^n[k])}. (8) 
(i,fc)ez 2 

Here the synthesis wavelet i/)(x; Tj[fc]) is derived from the mother wavelet ip(x) 
through the action the fHT corresponding to the shift Ti[k) — (f>i{k\/ir. The uni- 
tary nature of the fHT ensures that these fractionally-shifted synthesis wavelets 
have identical norms. In particular, while the amplitude \ci[k]\ indicates the 
strength of wavelet correlation, the local signal displacement gets encoded in 
the shift n[k] which specifies the most "appropriate" wavelet within the family 

2.2 Characterization of the Gabor-like transform 

It turns out that the shifted wavelets ip{x; Ti[k]) in (8) can be explicitly charac- 
terized when ip(x) is a (real) Gabor-wavelet, 

ip(x) = <p(x) cos (uj x + £ )- (9) 

This formula describes the asymptotic form of the Gabor-like wavelet ^(x; a) 
in (6). For reasons that will be evident shortly, we choose to flip the roles of 
the analysis and synthesis wavelets: we will analyze the signal using the dual 
complex wavelet W(x; a) = ip(x;a) + jtp'(x;a), while the Gabor-like wavelet 
^(x; a) will be used for reconstruction. 
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If the window function cp(x) is bandlimitcd to [—0,0] with O < ujq, we can 
make precise statements on the dual-tree representation in (8). To do so, we 
will need the following result: 

Proposition 2.1 Let <p(x) in (9) be bandlimited to (— Uo,u>o)- Then 

^t{(P(x) cos(ll> x)} = <p(x) cos(u> x + 7tt). (10) 

That is, the fHT acts on the phase of the modulating sinusoid while preserving 
the Gaussian envelope. In particular, we can then rewrite (8) as 

fixed window variable amp— phase oscillation 

f( x )= X! VtA x ) | c A k \ | cos (w a; + £o + [k]) | (11) 

(i,k)ez 2 

where <p% t k{x) — S^fecp(x) denotes the (fixed) Gaussian-like window at scale 
i and translation k. This provides an explicit interpretation of the parame- 
ter Ti[k] as the phase-shift applied to the modulating sinusoid of the wavelet. 
In effect, the oscillation is shifted to best fit the underlying signal singular- 
ities/transitions while the localization window (pi.k(x) is kept fixed. In this 
light, one can interpret the associated dual-tree analysis as a multircsolution 
form of the windowed- Fourier analysis, with the fundamental difference that, 
instead of analyzing the signal at different frequencies, it resolves the signal 
over different scales (or resolutions). 

Figure 1 shows quadrature pairs (^-^(x; a) , ^£ r+ i/ 2 ip(x; a)) of Gabor-like 
spline wavelets corresponding to different r. Each of the pairs are localized 
within a common Gaussian-like window, and the modulating oscillations are 
driven to a relative quadrature through the action of the pair (^-,^-+1/2)- 



3 BIVARIATE EXTENSION 

The amplitude-phase representation derived in §2 can also be extended to the 
2D setting where the dual-tree wavelets exhibit better directional selectivity 
than the conventional tensor-product (separable) wavelets [4]. 



3.1 Directional HT pairs of wavelets 

We briefly recall the construction framework for the bivariate DT-CWT based 
on the tensor-products of one-dimensional analytic wavelets [3] . Specifically, let 
<p(x) and cp'(x) denote the scaling functions associated with the analytic wavelet 
ip a (x) = ^(x) + jtp'(x), where tp'(x) = Jftp(x). The 2D dual-tree construction 
then hinges on the identification of four separable multiresolutions of L 2 (R 2 ) 
that arc naturally associated with the two scaling functions: the approximation 
subspaces V(y) ® V((p), V(<p) <g> V((p'),V((p') ® V((p) and V((p') ® V ((?'), and 
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their multiscale counterparts. The corresponding separable wavelets - the 'low- 
high', 'high- low' and 'high- high' wavelets - are specified by: 





= v{x)tp{y), 


ip4,(x) 


= <p{x)ip'(y), 




= i>{x)<p{y), 


$ 5 (x) 


= i>{x)<p'{y), 




= i>(x)4>(y), 


■tpe(x) 


= ip{x)i)'{y), 


■0 7 (x) 


= <p'{x)i>(y), 


4>io{x) 


= cpW(y) 




= ip'{x)<p(y), 


-ipu(x) 


= i>'(x)<p'(y) 


4> 9 (x) 


= ip'{x)ip(y), 







(12) 

The dual wavelets ^\{x), . . . , i/> 12 (a;) are similarly defined in terms of ip(x) & n d 
ip'(x) (here x = (x,y) denotes the planar coordinates). As far as the identifica- 
tion of the complex wavelets is concerned, the main issue is the poor directional 
selectivity of the 'high-high' wavelets along the diagonal directions. This prob- 
lem can, however, be mitigated by appropriately exploiting the one-sided spec- 
trum of the analytic wavelet ip a (x), and, in effect, by appropriately combining 
the wavelets in (12). In particular, the complex wavelets specified by 

#i(x) = ip a (x)(p(y) = i) 2 {x) + jtps{x), 
^ 2 (x) = ip a (x)(p'(y) = 4> b (x) + ji>u{x), 
* 3 (x) = cp {x)ip a (y) = tp^x) + j^i(x), 
^i(x) = (p'(x)ip a (y) = Vv(x) + jtpw{x), 

*,(*) = ±M*)Mv) - ( ^ x) ~^ x) ) +j (^ {x) + 



*.(») = ^ a {x)Mv) = { T2 T2 J, (13) 

exhibit the desired directional selectivity along the primal orientations 0\ = 
O2 = 0, 3 = 4 = tt/2, 5 = tt/4, and 6 = 3tt/4, respectively [3]. The dual 
complex wavelets ^>i(x), . . . ,^q(x) are specified in an identical fashion using 

the dual wavelets 'tpp(x), and are oriented along the same set of directions. 

Akin to the HT correspondence, the complex wavelet components are related 
through the directional HT (dHT): 

^ef(x) -jsign(u2 , w)/(w) (0 < < tt), (14) 

where ug = (cos 0, sin 0) denotes the unit vector along the direction 0. In 
particular, we have the correspondences 

3m(*/) = ^5Re(*/) (* = 1,...,6), 

so that, by denoting the real component of the complex wavelet ^fe(x) by ipe(x), 
we have the convenient representation ^e(x) = ipi(x) + jjfe e tp£(x) that is rem- 
iniscent of the ID analytic representation. 
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Figure 1: Quadrature pairs of Gabor-like spline wavelets obtained by the action 
of fHT group. Blue (solid line): Mrip(x;8), Red (broken line): H r+ iip(x;S), 
and Black (solid line): Common localization window given by \^Crifj(x;8) + 
jy£ r+ iip(x; 8)| (http:/ /dx.doi.org/doi. number. goes. here). 

3.2 Directional amplitude-phase representation 

Let us denote the dilated-translated copies of the each of the six analysis wavelets 
^(cc) by ^i,i,h(x), so that 

^i,i,h(x) = Ei^e(x) (ieZ,fce Z 2 ), 

where S^fe is specified by H$ k/(x) — 2 z f(2 l x — k). The corresponding dual-tree 
transform involves the analysis of a finite-energy signal f(x) in terms of the 
sequence of projections 

<t[k] = \(f,*l,i,h). 

The representation of f(x) in terms of the analysis coefficients cf [k] is based on 
the following fractional extension of the directional HT operator: 

Mq t = cos(7tt) J - sin(7rr) ffle (r e R). 

These operators allow us to capture the notion of direction-selective phase- 
shifts. The key properties of the fHT, which played a decisive role in establishing 
the representation for the ID counterpart, carry over directly to the fractional 
directional HT (fdHT) operators: they are invariant to translations and dilations 
and are unitary. In particular, based on the above properties, we derive the 
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representation 

/(*)= Y, \4[k]\^ k {M^r![k})} (15) 

(e,i,k) 

involving the superposition of direction-selective synthesis wavelets affected with 
appropriate phase-shifts [5]. The wavelets ipe(x; r/[fe]) are derived from the ref- 
erence wavelet ipe(x) through the action of fdHT, corresponding to the direction 
Qt and shift r/[fe] = arg(cf [fc])/7r As in the ID setting, further insight into the 
above representation is obtained by considering wavelets resembling windowed 
plane waves. 

3.3 Directional Gabor-like analysis 

Akin to the ID setting, further insight into the above representation is obtained 
by considering wavelets resembling windowed plane waves. A distinctive feature 
of the dHT, that comes as a direct consequence of (14), is its phase-shift action 
in relation to plane-waves: it transforms the directional cosine cos(w^:c) into 
the directional sine sm{ugx). Moreover, what turns out to be even more crucial 
in the current context, is that the above action is preserved for windowed plane 
waves of the form 

(p (x) cos( flujx). (16) 

In particular, as a straightforward directional extension of (10), we have the 
following generalization for the fractional extensions: 

Proposition 3.1 Suppose that <p(x) in (16) is bandlimited to the disk {u> : 
||w|| < fi}. Then we have that 

tig.T |cp(a;) cos(Qugx)} = <p(x) sin(flujx + ttt). (17) 

Thus, the fdHT acts only on the phase of the oscillation while the window 
remains fixed. In particular, if the dual-tree wavelets are of the form ipi( x ) = 
(pe(x)cos (fl(uj e x), we can then rewrite (15) as 

fixed window variable amp— phase directional wave 

' s f ' " * -I 

f(x) = Vt,i,k{ X ) S ^fe{ \4[ k }\ C0S ( n eUg e X + TTT^[k}) j, (18) 

(C,i,k) 

where (pe,i,k( x ) represent the dilated-translated copies of (pe(x). The above rep- 
resentation explicitly highlights the role of t\ [k] as a "scale-dependent" measure 
of the local signal displacements along certain preferential directions. This is 
the scenario for the spline-based Gabor-like transforms [3] where the dual-tree 
wavelets asymptotically converge to directional Gabor functions. 

4 CONCLUDING REMARKS 

We presented an amplitude-phase representation of the dual-tree transform in 
general, and a windowed-Fouicr-like characterization of the Gabor-like trans- 
forms in particular. The signal representation was centered around one crucial 
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construction, namely the HT correspondence between the wavelet bases. In- 
deed, the identification of the fHT-transformed wavelets in (11) followed as a 
direct consequence of this particular relation; the subsequent developments were 
then based on two crucial properties of the fHT, namely 

• its intrinsic invariances with respect to translations, dilations and norm- 
evaluations, and 

• its particular phase-shifting action on the Gabor wavelet. 

These observations could be of potential interest in applications involving the 
dual-tree transform, particularly signal denoising, where a rigorous mathemat- 
ical model linking the reconstructed signal to the processed complex wavelet 
coefficients is desirable. 

Appendix: Proof of Proposition 2.1 

Proof. The result follows directly from definition (7) and the following action 
of the HT: 

,¥C {(p(x) cos(u;ox)} = <P{x) sin(wox). (19) 

Indeed, we see that 

Hr{<p(x) cos(w a ; )} = cos(7tt) (p(x) cos(oj x) — sin(7rr) {(p(x) cos(cu x)} 
= cos(ttt) <p(x) cos(wox) — sin(7rr) (p(x) sin(uJox) 
= cp(x) cos(ojqx + 7tt). 

To establish (19), we note that the Fourier transform 1 of cp(x) cos(ujox) is given 
by 7r((p(w — uoq) + (p(ui + ujq))- Following definition (1), we then have that 

.3^ {(p(x) cos(u>ox)} < — ► — j sign(w) • 7r(cp(w — u>o) + (f»(w + wo)) 

= -jn(<(>{u - wo) - <P(^ + w )) 
<- — ► <p(x) sm(u) x), 

since — jir((p(ui — ujq) — cp(w + luq)) is the Fourier transform of (p(x)sm(oj x). 
Note that, in going from the first to the second step, we have used the crucial 
fact that the supports of (p(w + coq) and (p(o> — u>o) are entirely restricted to the 
half- lines {uj < 0} and {ui > 0}, respectively. 
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